Buffer gas induced collision shift for the 88 Sr 1 So — 3 Pi clock transition 
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Precision saturation spectroscopy of the 88 Sr 1 So — 3 -Pi is performed in a vapor cell filled with 
various rare gas including He, Ne, Ar, and Xe. By continuously calibrating the absolute frequency 
of the probe laser, buffer gas induced collision shifts of ~kHz are detected with gas pressure of 
1-20 mTorr. Helium gave the largest fractional shift of 1.6 x 10 - Torr -1 . Comparing with a 
simple impact calculation and a Doppler-limited experiment of Holtgrave and Wolf [Phys. Rev. A 
72, 012711 (2005)], our results show larger broadening and smaller shifting coefficient, indicating 
effective atomic loss due to velocity changing collisions. The applicability of the result to the 
So — 3 Po optical lattice clock transition is also d' 1 
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The latest progress of optical clocks is so rapid that 
their capability is surpassing that of microwave clocks. 
A fractional accuracy of ~ 10 -17 was demonstrated with 
single- ion systems [l| . This remarkable precision, on the 
other hand, requires us to consider various systematic 
shifts which have been so far ignored. Neutral atomic 
clocks, particularly fountain clocks, have been often lim- 
ited in their accuracy by binary collisions (2j. Besides 
that, whether the clock is neutral- atom or ionic system, 
collision shifts induced by background gas (BG) Q may 
ultimately limit the accuracy in this unprecedented level. 
However, the characterization of BG induced collision 
shifts is difficult in real clocks because the deliberate 
injection of buffer gas for the shifting of atomic reso- 
nance normally causes loss of atoms, limiting the coher- 
ence time needed to evaluate subtle BG induced shifts. 

From the beginning of laser spectroscopy, frequency 
shifts and broadening caused by background gas colli- 
sions have been investigated both experimentally and 
theoretically The spin-forbidden transitions 1 5'o— 3 Pi 
of atoms with two outer electrons have been major test- 
ing grounds due to their fairly narrow natural linewidths 
and level simplicity free from hyperfine structure. Ex- 
periments were done for Ca Q, Sr 0, 0, 0, Cd @, 1, 
Hg @, [Io| and Yb [U with perturber pressures of 
10° — 10 2 Torr. fn these experiments, shifts and broad- 
ening exceeds 10 MHz, and high order effects such as 
collision time asymmetry or speed-dependence of colli- 
sional broadening and shifting can be observed [H, [l3[ . 
The higher-order effects, however, disturb the accurate 
measurement of the linear collision shift and broadening, 
which is necessary to evaluate systematic shifts of atomic 
clocks. 

In this paper, buffer gas induced collision shifts for 
the 88 Sr 1 Sq — 3 Pi transition is reported based on pre- 
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FIG. 1: (Color online) Experimental configuration of the pre- 
cision saturated absorption spectroscopy. ECDL: extended 
cavity diode laser, FR: Faraday rotator, EOM: Electro-optic 
modulator, AOM: acousto-optic modulator 



cision saturation absorption spectroscopy with the FM 
sideband technique Q [l5j|. Continuously monitoring 
the probe laser frequency using a frequency comb refer- 
enced to International Atomic Time (TAI) , we compared 
the resonance frequency with the true atomic resonance 
measured with ultra-cold atoms [n| . The detected back- 
ground gas induced collision shift was ~ kHz, which is, to 
the best of our knowledge, the hnest definitive measure- 
ment of optical frequency shifts induced by foreign-gas 
collisions. 

The (5s 2 ) 1 S - (5s5p) 3 Pi transition at A = 689 nm 
has a natural line width 7 = 2tt x 7.4 kHz which corre- 
sponds to the saturation intensity / s = 2-7r 2 ?ic7/3A 3 = 
2.9 /LtW/cm , where c and 2-kK is the speed of light 
and the Planck constant, respectively. The experimental 
setup of our saturated absorption spectroscopy is shown 
in Fig. 1. The spectral width of an extended cavity diode 
laser (ECDL) is reduced to less than 500 Hz by lock- 
ing the laser frequency to a cavity with a basic Pound- 
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Drever-Hall (PDH) method [l7|. The modulation index 
and frequency u>m are 0.17 and 2n x 21.7MHz, respec- 
tively. This modulation is used for the saturation ab- 
sorption spectroscopy as well. Strontium is contained in 
a heat pipe with Brewster window plates at both ends. 
The typical temperature of the pipe is 728 K. The two 
ends of the active area (20 cm length) are water-cooled so 
that the strontium is not deposited on the windows. The 
buffer gas is introduced from one end of the pipe, and 
Pirani gauges are installed at both ends for the confir- 
mation of uniform pressures in the tube. The magnetic 
sublevels are well defined by a bias magnetic field of 3 
G produced by a Helmholtz coil pair. Pump and probe 
beam of the saturation spectroscopy have the same linear 
polarization which excites the 7r transition of Am,/ = 0. 
Pinholes are placed in both arms to overlap the pump 
and probe beam. The peak intensity and 1/e 2 radius of 
the pump beam (probe beam) are 270 //W/cm and 4.2 
mm (50 /jW/cm and 3.3 mm), respectively. After inter- 
acting with Sr vapor, the polarization of the probe beam 
is rotated 90 degrees and coupled out by a polarization 
beam splitter. The pump beam is frequency shifted by 
80 MHz from the probe beam in order to avoid DC noise 
caused by the interference between the two beams. In 
addition, the pump beam is chopped with a frequency 
of 5 kHz to extract the saturation signal induced by the 
pump beam. The signal is first demodulated at lom and 
then lock-in detected. The absolute frequency of the laser 
is always monitored by a frequency comb referenced to 
International Atomic Time (TAI). 

A sample of PDH phase detection signal with helium 
buffer gas of 4 x 10~ 4 Torr is shown in Fig. [UJa). The 
absolute frequency of the probe laser is expressed as a dif- 
ferential frequency (/ pump + / pro bc)/2 - /o, where / pump 
and /probe is pump and probe laser frequency, respec- 
tively, and the true atomic resonance /o is the frequency 
precisely measured in [l6j from ballistically expanding ul- 
tracold atoms. The full width half maximum (FWHM) 
is nonlinear- fitted to be 107 kHz, which agrees with satu- 
ration broadening (80 kHz) and transit time broadening 
(31 kHz) expected from the beam intensity and diameter. 
This dispersion shape works as an error signal to lock the 
laser. The stability of this Sr vapor cell clock is shown in 
Fig. [D(b) as modified Allan deviation against a Hydrogen 
maser. The bottom of the instability is 4 x 10~ 14 at 100 
s. The stability is limited by fluctuations of buffer gas 
and vapor pressure or pump beam intensity that changes 
the light-pressure- induced lineshape asymmetry [l8l |. 

We measured the pressure dependence of the back- 
ground gas induced broadening and shift for various 
background gases. The buffer gases are rare gases, 
namely helium, neon, argon, and xenon. The attempts 
to fill with nitrogen failed since the gas was quickly ad- 
sorbed by solid strontium. The broadening is determined 
from the dispersion curve as in Fig. [Ha). Locking the 
laser frequency to the zero-crossing of dispersion curve, 
the resonance is determined by a frequency comb that is 
referenced to TAI. We determined the slope of the colli- 
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FIG. 2: (a)Dispersion signal of FM side band spectroscopy. 
The offset frequency of f =434 829 121 312 kHz is the true 
atomic resonance measured in (b) Modified Allan stan- 
dard deviation of the 689 nm laser frequency tightly locked to 
the resonance signal obtained by the saturation spectroscopy. 



sional broadening and the frequency shift as a function 
of the buffer gas pressure by doing a least square fit, and 
the result is summarized in Table [T] 

The frequency shift data are summarized in Fig. [3] 
The experimental condition described above has resulted 
in a residual offset of -4 kHz from the true atomic reso- 
nance. Part of the reason is a shift caused by homonu- 
clear two body collisions. By changing the temperature 
of the heat pipe, we measured the Sr-Sr collision shift 
coefficient to be -1.3 x 10" 10 Hz cm 3 [lj|. This effect 
causes a shift of -1 kHz for the condition of Fig. 3. We 
suspect that the most probable reason for the remaining 
-3 kHz is residual Doppler shifts. The spatial filters in the 
two arms still permit a mutual beam tilt of 6 = 62 /xrad. 
In the case that atoms are directed as in an atomic beam, 
this tilt causes Doppler shift of ki,v m sin (6/2) = 15 kHz 
maximum, where ki is the wave number of the light, 
v m = y / 2ksT/m = 371 m/s is most probable speed of 
Sr atoms at gas temperature T = 455°C, m is the Sr mass 
and k,B is the Boltzmann constant. Our spectroscopy is 
ideally performed with an isotropic velocity distribution 
of atoms. However, a slightly biased velocity distribu- 
tion may be produced in the heat pipe due to an inho- 
mogeneous temperature distribution of the pipe. Since 
the beam passes above bulk solid strontium, the distri- 
bution could be biased upward. Other possibilities are 



light-pressure-induced lineshape asymmetry 18 1 and re- 
coils due to photon absorption and emission 2C , [2l|, [23] • 
Only the collisions with He show a positive shift coeffi- 
cient, whereas Ne, Ar, Xe show negative coefficients. The 
open symbols at lower pressures are excluded from the 
fitting since the vapor pressure of strontium is not small 
enough comparing with the buffer gas pressure. Even 
slightly above the vapor pressure, the buffer gas pres- 
sure wasn't uniform at pressures below 3 x 10 -3 Torr. 
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FIG. 3: (Color online) Shifts of the resonance frequency in- 
duced by various kinds of buffer gas. The corner of the slope 
around 3 mTorr is attributed to the inhomogeneous pressure 
distribution in the cell, which occurs due to insufficient buffer 
gas pressure relative to strontium vapor. Inset is the broad- 
ening of the Doppler-free atomic resonance. 

We added the buffer gas from one side of the pipe, and 
the spatial gradient of the gas pressure was observed by 
Pirani gauges placed at the both ends of the heat pipe. 

The collisional width T (HWHM) and shift A in the 
impact approximation are proportional to the density of 
perturbers N and are given by the well-known expression 

T+iA = N J d?v r f^v r )v r J dpp[l- {SaSj})^ } 

where / M (tv) is the Maxwellian distribution of relative 
velocities v r of colliding absorber and perturber hav- 
ing reduced mass /i, and p is the impact parameter. In 
this work we report results obtained in the mean speed 
approximation in which calculations arc done just for 
the mean relative speed v = ^SksT/ (irp) of colliding 
atoms. For the intercombination transition 1 Sq — 3 Pi, 
in the adiabatic approximation the angular average of 
scattering matrix elements in the initial and final states 

(SaSjfy = ±e-"> A - x + l e - ir > B - x can be expressed 

in terms of phase shifts r]A-x and T]b-x induced by the 
quasimolecular transitions A 3 + —X 1 + and B 3 \~ X 1 + , 
respectively [HI, H^j. These phase shifts are integrals 
?] = f_™ dtAV (r(t)) jh of the excited and ground elec- 
tronic state potentials differences V(A 3 + ) — V(X 1 + ) 
and V{B 3 l) - V(X 1 0+) for Sr-Rg(He, Ne, Ar, Xe) sys- 
tems. For simplicity we assumed straight line trajectories 
with the interatomic separation r(t) = \J p 2 + v 2 t 2 . 

The interaction potentials are approximated here by 
a modified Lcnnard-Jones formula V(r) = C i 2 /r 12 + 
C s /r s + C 6 /r 6 . Following Hindmarsh et al. [25[ C 12 
coefficients are estimated using Hindmarsh radii rjj (Rg) 



for rare gases Rg [2(| and rjj(Sr) for Sr. The Hind- 
marsh radii rjj (Sr) have been estimated usin g th e charge 
densities calculated by Mitroy and Zhang [271]. They 
were found to be 9.4773 ao and 10.8801 ao in the 
ground 1 S'o and excited 3 Po,i,2 electronic state of Sr 
atom, respectively. The atomic unit of length is ao ~ 
0.05292 nm. The dispersion coefficients C§ and C$ 
for Sr-Rg were also calculated by Mitroy and Zhang 
[27l ]. The dispersion coefficient given for excited states 
3 S + and 3 n were converted to appropriate coefficients 
for states A 3 + and B 3 1 using the following asymp- 
totic relations: C 6 (A 3 0+) = C 6 ( 3 U) and C 6 (B 3 1) = 
±C 6 ( 3 £+) + iC 6 ( 3 n). The term C 8 /r 8 included in our 
calculations allows for a better description of the interac- 
tion anisotropy in the excited state of Sr. We note that 
the addition of the Cg /r 8 and C10 / r 10 terms to the inter- 
action model changes the calculated broadening about a 
factor of two. The change can be even bigger in the case 
of the shift and for Ne can lead to a change of its sign. 
We also verified that if the thermal average is neglected, 
the resulting error is less than 5%. The positive colli- 
sional shift observed for He is caused by the dominating 
influnce of the repulsive interaction C12/?* 12 on the colli- 
sional broadening and shifting. 

Table |T] shows the comparison of the pressure broad- 
ening and shifting coefficients measured in this work 
with those calculated in the impact approximation, Eq. 
(fTj), and experimental results obtained by other authors. 
There is only one full set of data for rare gas pressure 
broadening and shifting of strontium intercombination 
line reported by Holtgrave and Wolf . Results of Ref. 
Q were obtained in the Doppler-limited regime at a tem- 
perature of 660(47) K, where the magnitude of the broad- 
ening and shifts exceed 100 MHz. Note that collisional 
broadening cross section determined by Holtgave and 
Wolf agree very well with those obtained by Crane et 
al [i[ for Ar and and differ by 23% from those for Ne. In 
Ref. [|[ also the Doppler-limited spectroscopical method 
was used. It is clear that our results do not agree with 
those obtained using Doppler limited spectroscopy [1, [|| . 
On the other hand our results agree well with those de- 
termined by Tino at al. Q for Ar in the Doppler free 
regime, see Tabic fl] 

The comparison of calculated results seems to give a bit 
better agreement with results obtained in the Doppler- 
limited case. Clearly Doppler-free results are charac- 
terized by bigger collision broadening coefficients and 
smaller magnitude of collision-shifting coefficients than 
those obtained from the Doppler-limited case. This can 
be attributed to velocity changing collisions [28| which 
change the absorber velocity significantly enough to push 
away the absorber from the condition of Doppler-free res- 
onance. Therefore, such events of collisions which can 
cause large frequency shifts in the Doppler-limited regime 
would not contribute to the signal in the Doppler-free 
case. Being a loss mechanism, it will lead to broadening 
but not to shifting of the line. This effect needs fur- 
ther theoretical and experimental investigation. More- 
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TABLE I: The comparison of our experimental pressure broadening and shifting coefficients with those calculated in the 
impact approximation, as well as with experimental results obtained by Holtgrave and Wolf [f| in the Doppler limited regime 
and obtained by Tino et al. [f| in the Doppler free regime. All data are given in MHz/Torr and the one standard deviation 
uncertainties are given in parentheses. 
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r/p 
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He 


3.79(17) 


+1.23(6) 


3.89 


+0.68 


8.67(10) +0.68(8) 


3.68 


+0.66 


Ne 


2.60(13) 


-0.64(4) 


1.17 


-0.31 


3.76(5) -0.06(3) 


1.10 


-0.28 


Ar 


2.45(4) 


-1.45(4) 


3.30 


-1.54 


4.76(8) -0.22(2) 


3.05 


-1.47 
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Xe 


3.44(23) 


-1.51(3) 


3.35 


-1.36 


4.71(6) -0.12(1) 


3.17 


-1.33 



"Result obtained by Tino et al. 0]. 



over in these investigations it would be very useful to 
carry out proper close-coupling broadening and shifting 
calculations [29j based on ab initio potentials. 

The cold strontium system is currently studied exten- 
sively in many national standards institutes due to its 
application to optical lattice clocks [30]. For the clock 
transition 1 Sq — 3 Pa, the interaction of the ground state 
rare gas atom with the excited clock atom at 3 Pq state 
is described by just one potential curve denoted 3 _ ap- 
proaching dissociation limit. The dispersion coefficient 
describing this interaction fulfills the following relation 
Ce( 3 0~) = Ce( 3 H) and with a good approximation is 
equal to Cq(A 3 Q + ) investigated here. Moreover, in the 
case of collisions of clock atoms with background gas 
atoms, the velocity-changing collisions will also play a 
crucial role leading to relatively large collisional broad- 
ening and small collisional shifting. Therefore, results 



obtained here as well as the discrepancy of the shift co- 
efficients depending on Doppler-free or -limited system 
should help to correctly estimate the impact of residual 
gas induced collision shift on optical atomic clocks. The 
values of fractional shift coefficients < 2 x 10 _9 Torr _1 
obtained here can be treated as an upper limit for shifts 
caused by the residual gas in optical lattice clocks which 
in fact should be even smaller. 
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Laboratory FAMO in Toruh, Poland. 
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